Are the Tonks regimes in the continuum and on the lattice truly equivalent? 
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Motivated by recent experiments, we compare the Tonks (i.e. hard-core boson gas) regime achieved in an 
optical lattice with the Tonks regime of a one-dimensional Bose gas in the continuum. For the lattice gas, we 
compute the local (i.e. on-site) two-body correlations as a function of temperature and the filling of the lattice. It 
is found that this function saturates to a constant value with increasing temperature. Furthermore, the parameter 
that characterizes the long-distance correlations in the lattice Tonks regime is also obtained, showing that on the 
lattice the long-distance correlations enter the Tonks regime more rapidly than in the continuum. 
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Strongly interacting gases of bosons in one dimension have 
been recently realized using optical lattices ul|j,|3|]. By load- 
ing a Bose-Einstein condensate into a deep two-dimensional 
optical lattice, an array of one-dimensional atomic systems 
(tubes) was created yllj]. Strong correlations amongst the 
bosons were subsequently induced by turning on a third lattice 
along the axis of the tubes, which further decreases the ratio of 
kinetic to interaction energy |l],|^. These anisotropic optical 
lattices exhibit a rich phase diagram " ^- In particular, by 
making the third lattice deeper and reducing the filling of the 
lattice below one particle per site, the Tonks regime, where the 
bosons effectively become hard-core and behave in many re- 
spects like fermions, was reached in the experiments reported 
in Ref. 2. Time of flight measurements of the momentum dis- 
tribution showed good agreement with a fermionization ap- 
proach which accounted for finite temperature, finite-size, and 
trap effects |2]. 

The achievement of the Tonks regime in an optical lat- 
tice raises a number of questions about the equivalence of 
this system with the continuum Tonks regime. Mathemati- 
cally speaking, the two types of Tonks gases correspond to 
the strongly interacting limit of two physically different mod- 
els. The lattice Tonks gas (LTG) is obtained from the Bose- 
Hubbard model (m = 1, . . . , Mq), 
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in the regime where 7^ = U/J ^ 1 1 19] and the filling of 
the lattice /q = Nq/Mq < 1 (iVo being the number of atoms 
and AIq the number of lattice sites). However, the continuum 
Tonks gas (CTG) regime is obtained from the Lieb-Liniger 
model M, 
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when the parameter 7 = Mg/fi^po ^ 1 |6, 7]. The pres- 
ence of the density po = Nq/L in this parameter is a dis- 
tinct feature of the continuum model. It means that the CTG 
regime can be reached either by increasing the interaction cou- 
pling g or by decreasing the density, po- By contrast, the LTG 
regime can be reached by increasing the ratio ji = U/J at 
any value of the filling /q. For low temperatures and fillings 
(i.e. /o ^ 1) the two Tonks regimes coincide since the Lieb- 
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FIG. 1: Local two-body correlation function gkiT^fo) = 
./"o"^((^m)^(fem)^> (times the square of 7l = U/J > 1 ) in the 
lattice Tonks regime, as a function of temperature T in units of 
the hopping J, for several lattice fillings /o < 1. At temperatures 
T ^ J, (/2^(r, /o) saturates to the value at half-filling, which is 
temperature independent. This behavior should be contrasted with 
the monotonic increase at high temperatures of the corresponding 
function for the continuum (i.e. Lieb-Liniger) model: giiT) = 
Po^{{-^\x)f{'i!{x)f) oc T (see Refs.HllJ. The dashed lines 
correspond to the low-T analytical approximation, Eq. <15> . 



Liniger model emerges as a low-density limit of the Bose- 
Hubbard model (see e.g. Ref. |8]). However, in |8] it was 
shown that, provided one has at most one particle per site (i.e. 
n„i < 1) and U ^ J, the Bose-Hubbard model can be effec- 
tively replaced by the following interacting fermion model: 

i^F = -^^(4 + iC„,+H.C.)+i/i+i/2, (3) 
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The couplings Ai = A2 = 7^ — J/U are small, and the 
fermions are 'almost' non-interacting. In this sense, therefore, 
one can speak of 'fermionization' of bosons. In this paper, we 
consider the regime where 7 or 7^ are large, which should be 
of interest for current 1 1, 2, 3] and future experiments explor- 
ing these correlated systems. 

The physical differences between the Tonks regimes of 
the Lieb-Liniger and Bose-Hubbard model can only be ad- 
dressed by an explicit calculation of their correlation proper- 
ties. These fall into two classes: short distance correlations 
are non-universal {i.e. model dependent), and therefore are 
expected to be different in the LTG and the CTG regimes. 
On the other hand, long-distance (or small momentum) cor- 
relations are characterized by the same power-laws exactly 
at the Tonks limit: 7 or 7^ -^ +00. Thus, for instance, 
in the thermodynamic limit at zero temperature, the momen- 
tum distribution n{p) ~ p^^/^ for p <$^ po 1^ [lOll . Nev- 
ertheless, for finite values of 7 or 7^, long-distance correla- 
tions are characterized by the Luttinger-liquid parameter K 
(e.g. n{p) ^ p^/^^~^), which is non-universal II2I1 . Pre- 
cisely, these non-universal features in the two Tonks regimes 
are what interests us here. In this regard, it is important to 
notice that for the Bose-Hubbard model there is no exact [i.e. 
Bethe-ansatz) solution available 1 13], and therefore analytical 
results for non-universal properties are scarce. In what fol- 
lows, we have summarized our results: 

• In the LTG regime we have obtained the temperature 
and filling-fraction dependence of the on-site two-body cor- 
relation function 33^ (T, /o) = fa^ {{bD^ibmf)- This 
is the lattice counterpart of the continuum-model g2{T) = 
Po^((vE't(x))2(*(a;))2). Both functions vanish for 7, 71, -^ 
+00, that is, when fermionization is complete. Whereas 
(72 (r) was computed in Ref. 1 14| for all 7 values, to the best of 
our knowledge no results existed for g2{T, /o). For 7l ^ 1 
and temperatures T -^ U,we find: 
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where f2{T) = {cm+2'^m) (^^^ below for details). At half- 
filling, the non-interacting fermion system is invariant under 
particle-hole symmetry: Cm -^ (~l)'"cj„^ which implies that 
/2(r) = at all temperatures. Thus we obtain the result 
that .g|'(T, /o = 1/2) = 27^^, independent of T. Results 
for arbitrary temperatures and several fillings are shown in 
fig-ID Although §2 {T, /o) is not directly related to the photo- 
association (PA) rate in the lattice because of the overlap be- 
tween the Wannier orbitals at different sites, it should be ex- 
perimentally accessible by suddenly ramping up the optical 
lattice before PA is performed in a time scale shorter than the 
atom tunneling time. After a substantial increase of lattice 
depth, overlap between the Wannier orbitals should become 
negligible. 

• We have obtained the Luttinger-liquid parameters K and 
Vs for fillings /o < 1 to leading order in 7^^: 

K ~ 1 + 47-1 sin ^/o/tt, (7) 

Vs/VF ^ 1-47^1 (/o COS tt/o). (8) 



where the Fermi velocity vp = Ja sin tt/q/S. Thus we con- 
clude that the Tonks regime is more easily reached, as far 
as long distance correlations are concerned, on the lattice 
than in the continuum. To see this, consider for instance a 
half-filled lattice (/o = 1/2). Using the above formula for 
7i ~ 10, K ^ 1.13, whereas for the Lieb-Liniger model 1121 
iiT ~ 1 + 4/7 = 1.4 for 7 = 7l = 10 (Indeed, 7 = 10 seems 
harder to achieve experimentally than 7^ = 10). The fact that 
long-distance correlations rapidly become Tonks-like for rel- 
atively shallow lattices justifies the fermionization treatment 
used in Ref. |2]. 

Next, we provide further details on the derivation of the 
above results. The key point is to notice that for temperatures 
T <^ U and filling less than one particle per site, the Bose- 
Hubbard model, Eq. Q, can be effectively replaced by the 
fermionic model of Eq. Q, ^f- In particular, the replacement 
can be made for computing the (low-temperature) partition 
function of Q: 
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where [i = l/T and iJjnt = Hi + H2 (cf. Eqs. g] 
|5})- In the second expression Zo = Tre^'^^^"^^^) and 
(...) = Tr[po{M,/9)---], withpo(Ai,T) = e-P^"^-^""^ /Z^, 
and Hq — Hy — Hint ■ We can now expand Z to the the low- 
est order in 7^^ and obtain |20], log(Z/Zo) = -/9(i?int) + 
0(7^^). Thus we need to compute the thermal average of 
i/jnt, which can be readily done with the help of Wick's theo- 
rem. The result can be written as follows: 
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where we have denoted {I ~ 0, ±1, ±2): 
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the function n(e, z) = \z' 



-l^pe 



1] is the Fermi-Dirac 



distribution for a fermion gas of fugacity z 
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— Jcos(pa) is the single-particle dispersion. Interestingly, all 
the above results follow from this simple expression, Eq. ( I10> . 
We begin by describing the calculation of q^ {T, /o). In an 
analogous manner to the continuum case 1 8, 14], this function 
can be obtained using the Hellmann-Feynman theorem: 
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By setting Ai = A2 = 7^^ in ( I10> . and assuming periodic 
boundary conditions ]21] so that /; = /_;, one obtains the 
first result given above, Eq. (|6}. An alternative expression for 
g^(r, /o) can be obtained after recasting 



/2(T) = /o(2J-2e2(T,^)_l 



(14) 



wher e e^(r, z ) = f„'^ Jde e'^g{e)n{e, z), being 5(e) = 
I/ttV"/^ — e^ the single-particle density of states. The advan- 
tage of this form of /2 (T) is that the Sommerfeld expansion 
can be used to extract the low-temperature behavior: 



52'^m/o)=5f(T = 0,/o) + 
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and52^(T = 0,/o) = 27^^ (1 - sin(27r/o)/27r/o). It is 
worth noticing that in the the low-filling limit /o -^ one re- 
covers, from the above expression, the asymptotic expression 
for the Lieb-Liniger gas obtained in Ref. [ij (see also jal), 
provided one makes the following identifications between the 
parameters of both models, M -^ f? j Ja^ {i.e. the effective 
mass), g — > Ua, po — > fo/a (see ||8|]). 

Finally, let us discuss some interesting properties of 
fo{T,z) and /2(T, z) defined in ( lll> . and their implica- 
tions for g2{T, fo). The first property is a consequence of 
the particle-hole symmetry of the non-interacting spectrum, 
which imphes that /o(T, z) + fo{T,z^^) — 1, and hence 
that the fugacity for filling 1 — /o is the inverse of the fu- 
gacity for filling /q. Likewise, one can show that f2{T, z) + 
/2(T, z^^) = 0, which implies that in practice it suffices to 
compute f2{T, z) for fillings /o < 1/2. The other property of 
these functions explains the saturation of .g^(r, fo) with in- 
creasing temperature observed in fig.^ for T ^ J, it can be 
shown that z -^ /o/(l^/o) and /2(T, z) -^ 0; hence the local 
two-body correlation function 32 (r, /o) — > g2{T, fo = 1/2) 
at all fo < 1. In the end, this is a consequence of the finite 
number of degrees of freedom available on the lattice. 

We finally consider the non-universal aspects of long- 
distance correlations, which are parametrized by the 
Luttinger- liquid parameters K and Vg . In order to obtain the 
latter, it is convenient to work with the density and phase stiff- 
ness iiir 
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Eq. ( I10> it can be seen that the terms responsible for the vio- 
lation, that is, for the renormalization of vj away from vp are 
those those coming from Hi . Thus, the renormalization of vj 
becomes manifest after noticing that [Hi , n„i] ^ 0, and there- 
fore one expects a non-zero contribution to the coefficient of 
dxj{x, t) (j{x, t) being the long wave-length part of the cur- 
rent density) in the coarse-grained continuity equation flSl]. 
This is one notable feature of Hp, which is obtained by pro- 
jecting on a low-energy subspace where n^ < 1 j^, in what 
may be regarded as a first step of the renormalization group. 

To sum up, we have shown that the Tonks regimes in the 
continuum and on the lattice are not, strictly speaking, equiv- 
alent. The local two-body correlations of the system on a lat- 
tice saturate with increasing temperature while for the Tonks 
regime of the continuum model are known to increase mono- 
tonically |8, 14]. Furthermore, the parameter K, characteriz- 
ing the decay of long-distance correlations Ilall51 . is more 
easily tuned to the Tonks limit on the lattice than in the con- 
tinuum. Finally, our results can also be extended to the cal- 
culation of corrections to the internal energy and entropy of 
the lattice gas 1 18]. However, the distinct behavior of the two 
Tonks limits is well displayed by the properties considered in 
this work, and we shall not pursue this task here. This research 
has been supported through a Gipuzkoa fellowship granted by 
Gipuzkoako Foru Aldundia (Basque Country). 

• Appendix: On the equivalence of the first and second 
quantization approaches to fermionization. The connec- 
tion between fermionization in the wave function formalism 
and its second quantization version, which leads to effective 
Hamiltonians like Eq. (|3}, has not been sufficiently empha- 
sized in previous treatments {e.g. Ref. 8). For completeness, 
we include a proof of their equivalence in this appendix. 

In their pioneering 1928 paper on second quantization of 
fermion fields, Jordan and Wigner flo] introduced a transfor- 
mation from hard-core bosons (or Pauli matrices) to fermions: 
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The angle <p corresponds to a twist in the boundary conditions: 
Cm+Mo = e'^c,„. This makes f^i{T = 0) ^ /,(T = 0) 
for / 7^ 0, and an explicit evaluation at finite size and zero 
temperature yields: 



/Kr = o,<^) = 
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Using that K = \/vj /vjm and Vg = y/VNVj lTj,[i3 and, to 
leading order in 7^ \ Eq = {G\Ho\G) + (G|Fi„t|G), where 
I G) is the ground state of Ho, along with ( II 6> . in the Mq -^ 00 
Eqs. and ^ are obtained (the same expressions were also 
derived by carefully taking the field-theoretic continuum limit 
ofHF 118]). 

Before concluding, an important difference between the 
Bose-Hubbard and Lieb-Liniger models is worth discussing: 
whereas the latter displays Galilean invariance, which implies 
that vj must be equal to the Fermi velocity vp ] 12], in the for- 
mer this symmetry is broken by the lattice. By inspection of 
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The operator Km (often referred to as the Jordan- Wigner 
string) turns the (hard-core) boson operator b„i into the 
fermionic c„j by attaching to it a phase factor which is de- 
termined by the number of particles to the left of site m. The 
trick converts the commutation relations of the 6's at differ- 
ent sites into the anti-commutation relations of the c's. In this 
appendix, it is shown that the same trick yields the celebrated 
Bose-Fermi mapping due to Girardeau IITII . 
the N-particle bra: 



Let us consider 
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where |0) is the empty state. The wave function 
$i?(xi, . . . ,xn) is anti-symmetric under exchange of any 



pair of coordinates as a result of the anti-commutation of 
the c's. This may lead us to think that the above bra de- 
scribes a system of N fermions. However, by noticing that 
K^^ ~ K„i, we can invert ( I18> and write the product 
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Next we shift all the string operators to the right and use that 
Krnfi) = |0) (since n;|0) = for all every time a string 
operator hits the empty state. Nevertheless, when commuting 
a string operator with a creation operator one must take care of 
a phase factor: KJ^^ = ^t-^e(x^-x^) y^^j^^ (where 61(0) = 
is assumed). After shifting all the string operators to the right, 
a factor like this one appears for each pair of particles, and 
therefore, 

where the fully antisymmetric prefactor A{x„i^ , ■ ■ ■ , Xm,^) = 



gi7ri:.<,e(x„. x^^) ^ n.<jSgn(a;m. - x^,)- Introducing 
the last expression into j20> . the bra can be rewritten as 
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^B y-^mi 1 ■ ■ • 1 ^rriN ) — -^['^mn • ■ ■ ■} ^mjsf )^F [•^mi ? ■ • ■ 5 •^m.N ) 

= |$F(a;mi,..-,a;™„)| (24) 

is a symmetric function which vanishes if x^m — Xm for 
i 7^ j. In other words, $b is the wave function of a system of 
hard-core bosons. This proves the equivalence of the first and 
second quantization approaches to fermionization. 
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